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The recent high-quality Boomerang data allow to test many competing cosmological models. Here 
I present a seven-parameter likelihood analysis of dark energy models with exponential potential and 
explicit coupling to dark matter. Such a model is conformally equivalent to a scalar field with non- 
minimal coupling to gravity. So far, the constraints on a dark energy - dark matter coupling were 
extremely weak. The Boomerang data constrain the dimensionless coupling /3 to be smaller than 
0.1, an order of magnitude better than previous limits. In terms of the constant £ of non-minimally 
coupled theories, this amounts to £ < 0.01. 
^ ' On the other hand, Boomerang has not enough sensitivity to put constraints on the potential 

slope. 
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I. INTRODUCTION 
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J> ■ In the last few years, a new reference model of structure formation has emerged, one in which 70% or so of the total 
matter content of the universe is in a form of dark energy. The main evidence for such a component is the supernovae 
la observations of Ref. jl| and J2|, which can be explained by assuming an accelerated expansion due either to a 
' cosmological constant or to a new matter component with equation of state 



p=(w-T)p (1) 



with w £ (0, 0.6) |jj H . This new component can also be modeled as a light scalar field with a potential that allows 
• for a potential-dominated epoch Q Q as, e.g., an exponential potential or an inverse power law . 

The evidence for the new component, sometimes denoted dark energy or "quintessence" Q, has been reinforced by 
the recent Boomerang CMB observations that show a preference for a flat universe (|] |t]] . In fact, since a CDM density 
of f2 c € (0.2, 0.4) is in agreement with a host of independent observations, from cluster abundance to cluster X-ray 
temperature, lensing, velocity fields etc., we can conclude that a conspicuous fraction of the total energy density has 
to be unclustered (or weakly clustered) and with negative pressure || , although the direct constraint on the amount 
J> , of dark energy from CMB alone is weak 

The nature of this extra component is so far completely unknown. The simplest explanation, the cosmological 
constant, runs against the argument that, from a dimensional point of view, its value should be more than a hundred 
orders of magnitude smaller than expected, and tuned with astonishing precision in order to become dominant just 
today. The hypothesis of some fundamental scalar field has at least the advantage that one may hope to build a 
theory that explains the coincidence as a consequence of some fundamental principle. For instance, the inclusion of 
a coupling between dark matter and dark energy, as will be done in this paper, might explain why the two energy 
densities are comparable. Another motivation for considering a scalar field lies in the fact that it is premature to 
be too specific about the dark energy properties: a general scalar field includes as a limiting case the cosmological 
constant, but allows also to investigate less extreme effective equations of state. 

In the same spirit of generality, I investigate the effect of an additional degree of freedom, represented by the 
coup ling of the scalar field to ordinary matter. Such a coupling has been proposed and investigated in refs. || |l(]] 



l[ [|2| 13 |l4| , and is equivalent, up to a conformal rescaling, to the classical Brans-Dicke coupling to gravity 



(for gravity coupling in the context of quintessence see refs. |T|] |T(|] |17| |L8| |L9|]). If the dark energy is coupled to 
dark matter alone (and not to baryons), the present constraints on the coupling are rather weak [^(J J2lj |l3|j ; as a 
consequence, in [^2| it was shown that the effects on the CMB turn out to be at the level of detectability already with 
the present data set. 

Aim of this paper is to use the most recent CMB data, the Boomerang power spectrum |jj , along with the COBE 
data (as reduced in (23)), to further constrain the coupled dark energy model. As particular cases, we will derive 
constraints on the pure cosmological constant model and on the uncoupled dark energy. A comparison of Boomerang 
with a different model of uncoupled dark energy is in [M. 

In all the calculations of this paper a flat universe is assumed. This is of course a severe limitation, but the number 
of free parameters is already large enough to be at the limits of computational capabilities. Moreover, the hypothesis 
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of quintessence has always been formulated in the context of flat models, in order to be consistent with the inflationary 
expectations. 

While the work was almost completed the results of the Maxima experiment [^6| pTj p8j have been published. 
They seem to confirm the results of Boomerang, but have not been included in the present paper. 



II. COUPLED DARK ENERGY 

The model of coupled quintessence or dark energy has already been studied in detail in |2^] . Here I limit myself 
to a summary of its properties. 

Consider three components, a scalar field <f> , baryons and CDM, described by the energy-momentum tensors 
2^1/(0) j T)ny(h) an d T)w(c)j respectively. General covariance requires the conservation of their sum, so that it is possible 
to consider a coupling such that 

T v(<l>);H = ( C cT( c ) + C b T {b) ) <f>. v , 

T u(b)^ = -C b T (b) <p- v . (2) 

where C c ^ b are the coupling constants. This particular coupling is indeed obtained by conformally transforming a 
non-minimally coupled gravity theory or, following [ p0[ , by the Lagrangian (in units G = c = 1) 

Ltot = + ~ U(4>) + L c (e 2C ^ 9lllJ ) + L b (e 2C ^g^) (3) 

where k 2 — 8tt and where L C:b denote the Lagrangian of the dark matter and baryonic fields. Notice that our constant 
C c corresponds to the invisible coupling constant /3j of [2(J. The radiation field (subscript 7) remains uncoupled, 
since T( 7 ) = 0. We derive the background equations in the flat conformal FRW metric ds 2 — a 2 (—dT 2 + Sijdx t dx : >) 
assuming the exponential potential 

U = Ae s<t > (4) 

as proposed e.g. in [p9| pi (3^]. As anticipated, we will couple the dark energy scalar field to the dark matter 
only, putting C b — 0. We call this choice dark-dark coupling. Generalizing |32| we introduce the following variables 
(putting H — a /a) 

K (f> na U na fpZ na fp b 

Adopting the e-folding time a = log a we can write the field, baryon and radiation conservation equation as a system 
in the four variables x, y, z, v that depends on the parameters (X, /? 

1 ) x — fxy ' 



1 )x - [iy + /3(1 - x 2 - y 2 - z 2 - v 2 ), 



fxxy + y\^2 

(l - 3x 2 + 3y 



v> = - V -(~3x 2 + 3y 2 ~z 2 ) (6) 

where the prime denotes derivation with respect to a, and where (3 = C c \J , /i = s^J ^y. The coupling constant 
j3 can therefore be regarded as the ratio of the dark energy coupling strength to the gravitational strength. The dark- 
dark coupling (3 is unconstrained by local experiments or by G/G measurements |2(]]. Constraints can be derived only 
by global effects on the cosmological dynamics: for instance, from the requirement of a universe older than 10 Gyr 
Ref. @ found 

1/91 < 1 (7) 
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In the following we show that Boomerang puts a much more stringent constraint on the dark-dark coupling. 

Neglecting the baryons, the system has a transient attractor (properly speaking, a saddle point) at the point 
xq = 2/3/3, uq — 0, zq = on which the scale factor expands as a ~ £ 6 /( 9 + 4 ^ ) (here t is the ordinary time defined as 
dr = a(t)dt), that is, slower than in a ordinary MDE. During this stage, first found in Ref. jjl), and denoted (/>MDE 
in p2| , the scalar field kinetic energy and the dark matter have a constant density ratio. When the universe leaves 
the saddle, it reaches a global attractor if \i < 3, at the point xo = —fi/3,yo = (1 — /i 2 /9) 1 / 2 , z — and the universe 
expansion follows the law a ~ t 3 ^ mimicking a perfect fluid with equation of state p = (w<j, — l)p, with — 2/i 2 /9 
. The expansion is accelerated only if /i < v3- As long as the baryon component is small with respect to the CDM, 
the phase space of the system remains qualitatively the same. 

It can be noticed that the 0MDE depends only on 0, while the subsequent (^-dominated epoch depends only on 
(X. Therefore, as shown in |33|, the possibility to set constraints on the coupling depends on the existence of the 
0MDE. This epoch is actually quite more general than appears from the discussion above: it exists infact for any 
potential such that U(4>) and dll/dcj) go asymptotically to zero, as for instance when U is an inverse power law. 

The coupled quintessence with exponential potential contains the case of pure cosmological constant (// = 0, = 0), 
of uncoupled quintessence {[3 = 0) and is asymptotically equivalent to a perfect fluid with a constant equation of state 
Wcf, — 2/i 2 /9. The model is conformally equivalent to a non-minimal gravity theory with Lagrangian term —^£ip 2 R 
and potential V(ip) = Xtp n with £ = 2/3 2 /3 and n = 4 + fj,/0 @ @ @. Therefore, a bound on amounts to 
constraining the non-minimal gravity coupling £. Notice that in the case we study here, in which the dark energy 
couples differently to baryons and to dark matter, there are two conformally related metrics in the Jordan frame 
(the frame in which gravity is coupled to the scalar field): the constant £ couples the scalar field to the Ricci scalar 
expressed in the metric in which the dark matter, rather than the baryons, follow the geodetics (see PC)] ). 

The perturbation equations, derived in p3[ ], are integrated by the use of a purposedly modified CMBFAST code 
p4j| . In addition to the scalar field, baryons, CDM and radiation, the code includes also massless neutrinos. I choose 
adiabatic initial conditions, as suggested by inflation. The initial conditions for the background equations are found 
for each set of parameters by trial and error so as to give today , fl c and Hq as requested. This procedure is so 
time consuming with respect to the code without coupled scalar field that the assumption of flat space and a further 
reduction in the parameter space explained below turned out to be necessary. 



III. LIKELIHOOD ANALYSIS 



Our theoretical model depends on two scalar field parameters and four cosmological parameters: 

0,fj,,n s ,h,Q.b,&c (8) 

The remaining input parameters requested by the CMBFAST code are set as follows: T cm f, = 2.72QK, Yh c = 
0.24, N v = 3.04, r c = 0. The latter quantity specifies the optical depth to Thomson scattering, and is a measure of 
reionization. In the analysis of JjJ this parameter was also included in the general likelihood and, in the flat case, was 
found to be compatible with zero at slightly more than la . Moreover, in ref. |?J it is found that fixing r = has only 
a minor effect on the other parameters. Therefore here, to further reduce the parameter space, I assume r c to vanish. 
Two other approximations with respect to @] have been necessary: first, I did not include the cross-correlation between 
bandpowers because it is not available. Second, an offset log-normal approximation to the band-power likelihood has 
been advocated by p3|| and adopted by 0, but the ir-quantity necessary for its evaluation is not available. Since the 
offset log-normal reduces to a log-normal in the limit of small noise I evaluated the log-normal likelihood 

- 2IogL(a i ) = T f^(^)-^)] 2 (9) 

where Zt = logCV, the subscripts t and d refer to the theoretical quantity and to the real data, Cg are the spectra 
binned over some interval of multipoles centered on ti, af are the experimental errors on Zg d, and the parameters 
are denoted collectively as aj. A 10% calibration error is added to the experimental errors [Q. I also evaluated for 
comparison a likelihood gaussian in the variables Ci, and found that the results do not change appreciably up to 2a 
from the peak. 

Among the parameters that refer to the scalar field, we have already shown in |22| that fi has a negligible effect on 
the background solution at z 3> 1, since the equivalence time does not depend on it and, although jj, sets the speed 
of the present accelerated expansion, this has only a minor effect on the perturbation spectrum at decoupling. This 
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will be confirmed below. Also, in 22 we have shown that the dynamics of the system is insensitive to the sign of 3, 
since both the 0MDE and the final accelerated epoch do not depend on it. We will consider only (3 > 0. 

In order to compare with the Boomerang analysis I assume uniform priors as in JtJ, with the parameters confined 
in the range 8 G (0,0.16), /iG (0,2.1), n s G (0.7,1.3), h G (0.45,0.9), f2& G (0.01,0.2), fl c G (0.1,0.9), with 
the further weak big-bang nucleosynthesis (BBN) constraint Vl^h 2 < 0.05 (and of course the condition f2 c + f2& < 1) . 
I evaluate the likelihood also using the most up-to-date BBN limit Pq] , assuming a gaussian prior for fl^h 2 with 
mean 0.019 and la error 0.002 (I refer to this as the strong BBN prior). The range of fi includes all the values for 
which there is acceleration at the present. A further parameter, the logarithm of the absolute normalization of the 
Ci spectrum, is integrated out analytically with uniform prior. The same age constraints (> 10 Gyr) used in is 
adopted here. Note that here we do not have the problem of near-degeneracy of parameters that is considered in ref. 
. In fact, the degeneracy that exists for those values of f7 c , f^, h, and 8 that give an identical angular-diameter 
distance to the decoupling surface is completely removed by the combined constraints of flatness, BBN and the allowed 
range of h and fi. We also checked that the limits on n s and 8 are broad enough to contain the bulk of likelihood. 

A grid of ~ 10, 000 multipole CMB spectra Ce is used as a database over which I interpolate to produce the 
likelihood function. Since both the COBE and the Boomerang data are in fact binned over intervals of multipoles, I 
average in the same bandpowers the theoretical spectra for a correct comparison. Three cases will be distinguished: 
pure A (fj, = 8 = 0); dark energy ( (3 = 0); coupled dark energy. In Fig. 1 I report examples of multipole spectra 
obtained varying fj, and (3 and fixing the other parameters: {n s , Qfe, £l c , h} — { 1, 0.05, 0.3, 0.7}. The quantity plotted 
is actually (£(£ + l)Ce/2n) 1 ^ 2 fiK as customary. The strong decrease in amplitude of the acoustic peaks as (3 gets larger 
than 0.1 depends on the fact that for these values the onset of the </>MDE occurs before the decoupling. During the 
0MDE the fluctuations smaller than the horizon grow less than during MDE, so that the fluctuations on the acoustic 
scales are depressed relative to the larger scales. Moreover, the addition of the integrated Sachs- Wolfe effect at small 
multipoles decreases the normalization of the intrinsic fluctuations at decoupling pl| |22l . 

In Fig. 2 ( panels a — d, solid lines) I plot the likelihood function for the parameters n s , h, f2(,, f2 c (marginalizing 
over all the other parameters) in the case of pure A and weak BBN prior. They are reasonably well in agreement with 
the analysis of (their case P10). The means and variances for the pure A model are 

n s = 0.96 ±0.06, h = 0.73 ± 0.09, Q b = 0.056 ± 0.016, ft c = 0.5±0.2, (10) 

while the best estimates (peaks) are n s — 0.975, h — 0.7, Slf, — 0.05, Q c = 0.3 in agreement with Q (notice that Ref. 
quotes the likelihood maximum for clear detections, e.g. for fl^ and n Sl and the mean in the other cases). For the 
case of uncoupled dark energy, I plot in Fig. 2 the likelihood functions in short dashed lines, plus in the panel e the 
new parameter [i. As expected, there is almost complete degeneracy in the direction of /z, and the results of pure A 
remain very similar. If the Boomerang errorbars were reduced to one third, the likelihood would begin to show some 
preference for higher values of u, as shown in the same panel e (dotted line). The problem of a high baryon content 
from Boomerang (see e.g. ref. jjTj) is not alleviated by this model of dark energy. 

The likelihood functions for the coupled dark energy case are plotted in Fig. 2 as long-dashed lines. Due to the 
degeneracy along /u, we can simplify the analysis fixing this parameter to any value in the relevant range; we put 
fi = 0.25. It turns out that the previous results are quite robust also with respect to the coupling, except for a shift 
of n s toward smaller values (n s — 0.88 ± 0.06). In the panel / I plot the coupling (3 : at the 96.8% c.l. we obtain 

101 < 0.1 (11) 

Assuming the strong BBN prior we obtain the likelihood functions shown as dot-dashed lines in the same Fig. 2 
(shown only for the coupled dark energy case for clarity; the other cases are similar). As expected, smaller fif, and 
smaller h than before are now acceptable. Also, fi c moves to smaller values, in order not to decrease excessively the 
first peak. This effect, along with the decrease of n s in panel 6, is observed also in Q when the strong BBN is imposed. 
The likelihood for the exponential slope /i (not shown) becomes quite flatter with the strong BBN. Interestingly, j3 
now peaks around 0.05, because this value raises somewhat the first peak, compensating for the low baryon content. 
The limit (|ll| ) remains valid. 



IV. CONCLUSIONS 



The most interesting conclusion that can be drawn from our analysis is that the coupling 8 between dark matter 
and dark energy has to be smaller than 0.1, an order of magnitude better than previously |2Cf |, independently of the 
BBN constraint. Although I derived this limit in the particular case of an exponential potential, it extends to a much 
larger class of theories since the 0MDE is dominated by the kinetic, rather than potential, energy of the scalar field, 
as will be shown in another paper, 
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In terms of the coupling constant £ of non-minimal gravity theories we get roughly £ < 0.01 at the 96.8% c.l. . 

The likelihood of the other cosmological parameters is robust with respect to the addition of the coupling. 

The potential parameter [i, which sets the effective equation of state of the dark energy at the present, is not well 
constrained by the present CMB, since its effect is very recent on the cosmological time scale. 

It is of course possible to add further constraints from large-scale structure, supernovae la and from other CMB 
experiments. In |22j a constraint similar to (|ll| ) was found from cluster abundance, for a particular choice of the 
other parameters. Future data, especially high-resolution, high-coverage experiments on CMB, have the potential to 
strenghten the limit by at least another order of magnitude. 
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FIG. 1. Top panel: spectra C\ = {1(1+ l)C e /2iv) 1/2 /iK of uncoupled dark energy for several values of /i (solid: fj, — 0; dots: 
fi — 0.7; short dashes: /i = 1.4; long dashes: /i = 2.1). Bottom panel: spectra C\ of coupled dark energy for several values of 
(3 (solid: /3 = 0; dots: (3 = 0.08; short dashes: (3 = 0.1; long dashes: (3 = 0.12) 
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FIG. 2. Marginalized likelihood functions. Solid lines: pure A model. Short-dashed lines: dark energy. Long-dashed 
lines: coupled dark energy. Dot-dashed lines: coupled dark energy imposing the strong BBN prior. In panel e the dotted line 
represents the likelihood for fx if the Boomerang errorbars were reduced to one third. 
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